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Figure 1. 



Figure 2. 



Figure 3. 



Figure 4. 
Figure 5. 



Figure 6. 



Ocean temperature distribution predicted using an embedded 
mixed layer model within an ocean circulation model. Cross- 
section is aligned along the track of the model storm which 
was located at 450 km at the initial time and has moved toward 
the right at 5 m s - ^ and is at 1315 km after 48 h. 

3 

Cumulative percentage of rank-ordered 3 h u # , u # and Q 

3 

values versus cumulative percentage of observations based on 
24-, 15- and 23-year samples (September to December) at OWS P, 
V and N, respectively. 

Daily-averaged values of (A) sea-surface temperature, (B) 
mixed layer depth, (C) insolation (Q s ) and surface heat flux 

3 

plus back radiation (Q ), and (D) u* at OWS P during 1963. 

3 

Dashed lines represent the long-term mean (24 year) values. 
Similar to Fig. 3 except at OWS V during 1959. 

Model-predicted maximum daily depths of the well-mixed layer 
(solid), surface temperature minus 0.2°C (top of shaded), 
surface temperature minus 1.0°C (top of cross-hatching) and 
surface temperature minus 2.5°C (bottom of cross-hatching) 
during 1959 at OWS P. 

Wind speed, total heat flux and solar heat flux at 32N, 175W 
during 1977 extracted from FNOC atmospheric model analyses and 
12 h predictions . 
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Figure 7. 



Figure 8. 



Figure 9. 



Figure 10. 



Figure 11. 



Figure 12. 



Location of TRANSPAC ship-of-opportunity bathythermograph 
observations during March and April 1977 (provided by W. White 
and R. Bernstein) . 

il -2 

Heat content (10 cal cm ) relative to the 200 m temperature 
calculated along 175°W from TRANSPAC analyses in September 
(solid) and November (dashed) 1976. Vertical bars indicate 
the cumulative surface heat flux between 15 September and 15 
November 1976. 

Initial (circles) temperature profiles at (a) 38°N, 165°W and 
(b) 32°N, 175°W from September 1976 TRANSPAC analysis and mean 
predicted values at 10 m intervals for months of October 
(triangle), November (horizontal dash) and December 1976 
(cross). Verification data from December 1976 TRANSPAC 
analysis are given at 0, 30, 60, 90, 120, 150 and 200 m 
(diamond) . 

Model-predicted layer (a) depth and (b) temperature during 
1976 along 38N. Initial spacing between the traces at each 
longitude are 100 m and 2°, respectively. 

Composite of model-predicted mixed layer temperature (upper) 
and depth (lower) for 6 longitudes as in Fig. 10, relative to 
spring transition date (zero day). 

Composite of atmospheric forcing parameters relative to spring 
transition date as in Fig. 11. Wind speed, total heat flux 
and solar heat flux are extracted from FN0C atmospheric model 
analyses and 12-h predictions. 
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Figure 13. Model-predicted mixed layer depth and temperature as in Fig. 
10 except along 155W during 1976. 

Figure 14. Model predictions as in Fig. 9 except at 30N, 175W and for 

March (circle), April (triangle), May (horizontal dash) and 

June (cross) 1977. Verification data are from June 1977 

TRANSPAC analysis. (a) Unmodified surface heat flux and (b) 

-2 -1 

reduction of upward heat flux by 10 cal cm h 
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1. Introduction 



One often hears that oceanography is several years behind meteorology. 
Many developments in physical oceanography parallel or build upon the 
earlier work in meteorology because of the similarity in the dynamics. It 
is useful to assess the stage of development of the oceanographic endeavor 
by comparison with the progress in meteorology. The example to be 
illustrated in this talk^ is the development of a capability for ocean 
prediction using numerical methods. 

The theme of this talk is that the next decade will see the development 
of a dynamic ocean prediction system comparable to that of the numerical 
weather prediction system developed during the 1960s. Development of an 
oceanic prediction system is a multi-faceted problem. Here, we pursue 
some aspects of the necessary theoretical background, representation of 
physical processes, observational-support systems and justification for an 
ocean prediction system. We limit our discussion to open-ocean regimes 
and thus avoid the especially complex circulations in coastal regions. In 
addition, we limit our choice of examples of research-in-progress to our 
group at the Naval Postgraduate School. Our intention is to use these 
examples (those of other groups could also be used) as illustrations of 
the progress toward an operational ocean prediction system. By analogy 
with the evolution of numerical weather prediction, it seems evident that 
the building blocks have been laid for development of a comparable ocean 
system. Justification for the system lies in the deployment of 

^This paper is based on a seminar presented at the Geophysical Fluid 
Dynamics Laboratory, Princeton, New Jersey, on 3 November 1979. 
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antisubmarine systems the U. S. Navy, conduct of fisheries management as 
guided by the Department of Commerce and perhaps also in climate research. 

It is not our purpose here to discuss the gaps in our knowledge that 
must be filled if an operational model is to be produced. Much research 
and development, many talented workers and considerable computer resources 
will be required to accomplish the task. As during the development of 
numerical weather prediction, there will be exciting and rewarding 
research opportunities as observationalists and modelers work toward a 
common goal. 

2. Theoretical background 

The basis for numerical weather prediction (NWP) was established in 
Princeton, NJ, at the Institute for Advanced Studies (e.g. Platzman, 

1979), but the development of operational numerical weather prediction 
models required years of additional research. Certainly the success of 
atmospheric general circulation models (GCM) at the Geophysical Fluid 
Dynamics Laboratory and other places was an important step in 
demonstrating the feasibility of NWP models. Such factors as numerical 
stability over long integration periods and the representation of physical 
processes in the GCM were necessary building blocks. 

Continuing this analogy with atmospheric models, the numerical 
simulations of the global ocean (see review by Haney, 1979), as well as 
the eddy scales, would seem to provide a similar demonstration of 
feasibility for ocean prediction. However, nearly all of these models 
have been driven with constant or slowly-varying forcing. The ocean mixed 
layer serves as a buffer zone between the atmospheric forcing and deep 
circulation. To resolve properly the range of time scales induced by 
realistic atmospheric forcing, we will require an ocean model that 
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includes both mixing and advective processes. One-dimensional mixing 
models simulate a major fraction of the upper ocean response to atmos- 
pheric storms (cf. Camp and Elsberry, 1978). Consequently, the recent 
advances in mixed layer theory (see review by Garwood, 1979) are an essen- 
tial component. However, there are other examples in which advective 
effects can not be neglected. 

One of the chief problems in coupling the mixed-layer models and the 
ocean circulation models is the difference in time and space scales. 
Whereas mixed layer models typically use time steps of 1 h and vertical 
increments of 1 to 3 m, the oceanic GCM may use time steps of many hours 
and may have layers 10 to 100 m thick. The approach at the Naval Post- 
graduate School has been to embed the bulk, turbulent closure model of 
Garwood (1977) within the vertical structure of the oceanic GCM of Haney 
(1980). An example (see Adamec , et al 1978, 1979) of the simulations 
which may be achieved with such a coupled model is shown in Fig. 1. Only 
the upper 600 m of the solution is displayed. A total of 25 layers rang- 
ing in thickness between 6 m near the surface and 100 m has been used over 
the 1000 m depth. A case of hurricane-ocean interaction is chosen because 
of the strong advective and mixing effects to be expected. A stress pat- 
tern typical of a hurricane has been moved from an initial point at 450 km 
to 1315 km after 48 h. As expected for a hurricane moving at 5 m s , a 
large amplitude oscillation is set up on the thermocline. Isotherms within 
the thermocline are first displaced downward as the storm center ap- 
proaches. Following the storm passage, there is a rapid upwelling, with 
the 17. 2° C isotherm being raised from around 200 m to about 105 m. Regions 
of large horizontal temperature gradients are produced in advance and 
behind the cold upwelled water. These thermocline oscillations would 
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continue for several inertial periods before being damped. There is a net 
temperature decrease near the surface due to the upward heat flux to the 
storm and the mechanical mixing induced by the strong winds. Even though 
the wind speed is decreasing after the storm center passes, cooling 
continues as the upwelling brings cold water nearer the surface where 
mixing is effective. 

The purpose of showing such an extreme example as hurricane-ocean 
interaction is to demonstrate that realistic advective and mixing effects 
can presently be simulated in a coupled ocean model. A relatively fine 
grid and detailed atmospheric forcing would be required to attempt a 
prediction with real data. Before developing such a complex model, it 
will be necessary to begin with simpler models. The "first-generation" 
model will probably be a mixed layer model only. The purpose of these 
models will be to predict ocean thermal structure changes in response to 
atmospheric forcing on time scales of days to weeks. In the following 
sections, we will describe the physical processes, the atmospheric forcing 
and the ocean observations necessary to run these first-generation models. 
3. Representation of physical processes 

One of the first prerequisites for atmospheric prediction was the 
proper representation of the processes involved in extratropical 
cyclogenesis. By the 1960's, there was ample theoretical background for 
demonstrating the length and time scales required for prediction of this 
phenomenon. It was also important to develop a representation of the 
internal energy redistribution and frictional effects acting on these time 
scales. We still have much to learn about the release of latent heat in 
clouds of different scales and the interaction with the atmosphere 
boundary layer, especially in the tropics. There is little doubt that 
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attempting to develop an ocean prediction model will also uncover some 
physical processes which will likely require additional years of research 
to model properly. 

The first-generation ocean mixed-layer models will take advantage of a 

separation of time scales. That is, the vertical mixing processes on time 

scales of a few days dominate the advective effects over much of the 

ocean. Thus the principal physical process to be represented is the 

redistribution of the density structure induced by vertical mixing 

processes. If the upper ocean is well mixed to a depth Z = -h, the 

vertical temperature flux within the layer will be a linear function of 

depth. Given the surface heat flux, the problem is to determine the 

vertical temperature flux at the base of the mixed layer. As shown by 

Niiler and Kraus (1977), this vertical temperature flux can be written in 

terms of an entrainment velocity w^ at the mixed layer base. 

One of the important factors influencing w^ is proportional to the 

2 -1 

third power of the atmospheric friction velocity (u* = xp ). The 

3 

distribution of u* values is very skewed, as shown in Fig. 2. In this 
diagram, the three-hourly observations at ocean weather ships have been 
rank-ordered and accumulated (see details in Elsberry and Camp, 1978; 
Elsberry and Raney, 1978). For example, 50 percent of the smallest values 

3 

of u# throughout the entire record contribute only 15 percent of the 

total accumulated value. The remaining 50 percent of the total rate of 

working on the upper ocean by the wind occurs during a few large wind 

speed events, lasting less than 15 percent of the time. These are 

associated with the passage of extratropical cyclones. This type of 

distribution occurs during both winter and summer seasons. It is 

3 

remarkable that the same curve fits the u* distribution for the three 
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ocean weather ship locations (P: 50°N, 145°W; N: 30°N, 140°W; V: 34°N, 
164°E) shown, especially considering the large differences in the total 

O 

u* . The distribution of the upward heat flux at the ocean surface (Q g ) 
is not as skewed, as it follows the u # (or wind speed) distribution. The 
Qa is a part of the buoyancy flux that also contributes to the entrainment 
velocity at the base of the layer, if the heat flux is upward. The con- 
vective overturning of parcels contributes to the turbulent energy avail- 
able for mixing at the base of the layer. 

3 

Daily values of u* , solar flux (Q ) and Q at OWS P are shown in 
* s a 

3 

Fig. 3. As indicated above, the u# trace is characterized by a 

background of small values with much larger values of brief duration at 

3-4 day intervals. There is also a synoptic period variability in Q and 

s 

Q . The oceanic response to this atmospheric forcing is shown in terms of 
a 

the sea-surface temperature and mixed-layer depth. From the beginning of 

September until about 10 October 1963, the temperatures were higher than 

the long-term mean. A major fraction of the seasonal decrease in 

sea-surface temperature then occurs over the next 10 to 15 days, and the 

temperatures remain below normal for the remainder of the year. The 

daily-averaged, mixed-layer depth trace is about 10 m shallower than the 

long-term mean prior to 10 October 1963. The layer deepens rapidly during 

3 

the period of large u* values, and remains deeper than normal throughout 
November 1963. 

We also consider the periods of light winds to be important in 

determining the oceanic thermal response to atmospheric forcing. An 

3 

example of the distributions of u,_ , Q g and at OWS V during January to 

3 

August is shown in Fig. 4. The large variability in u_ v relative to the 
long-term mean is again shown. Note that there is a period from around 15 
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*3 

March to 15 April during which the mean u # decreases significantly. The 
summer-time values are considerably smaller, although there continues to 
be synoptic-period variability. During the same period that the wind 
speeds are decreasing, the solar flux is increasing. A net downward heat 
flux tends to oppose deepening of the layer. If there is insufficient 
wind mixing to maintain the depth against the stabilizing influence of the 
surface heating, the layer will become shallower. This tends to occur 
throughout the year during the low wind speed periods between the passage 
of storms, as shown in both Figs. 3B and 4B. Likewise, the period of maxi- 
mum daytime heating will lead to a similar shallowing of the layer if the 
wind speed is not large enough to maintain layer depth. A very rapid 
transition in mixed-layer depth occurs around 1 April 1959. Prior to this 
time, the depth was around 130m. During a single diurnal period, the 
layer depth decreased to 40 m. Although the passage of subsequent storms 
temporarily increased the layer depth, the values did not approach the 
winter-time values. The effect on the sea-surface temperature is to trap 
the heat flux in a shallower layer, and thus increase the temperature. If 
the layer is very shallow, the rate of temperature increase can be quite 
large. Elsberry and Garwood (1978) have shown that warm and cold anoma- 
lies in the upper ocean throughout the summer can result from early and 
late transition dates, respectively. 

The important feature of these data sets is that the oceanic mixed 
layer depth and temperature do not smoothly change in time. Rather the 
evolution tends to respond directly to the frequency and intensity of 
storm events. Rapid cooling and deepening occur during high wind speed 
periods associated with storms, especially if there is upward heat flux. 
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Warming and shallowing tend to occur during light wind speed periods with 
downward heat flux. 

It is important to determine whether an oceanic mixed-layer model can 
predict these changes in the ocean thermal structure given the atmospheric 
forcing. An example (Garwood and Adamec, 1980) of a long-term simulation 
with the 3 h forcing at OWS P is shown in Fig. 5. The Garwood (1977) 
model was started from the observed temperature on January 1. This is the 
only ocean thermal structure observation used, although the observed 
sea-surface temperature has been used to calculate the surface heat 
fluxes. The simulated isothermal depths during winter indicates consider- 
able variability as in Fig. 3. A rapid transition to a shallow depth is 
predicted about day 100, after which the layer remains above 65 m. Forma- 
tion of the seasonal thermocline and its associated stability (temperature 
intervals relative to the surface temperature) are shown in Fig. 5. The 
autumn storms then gradually erode the layer. Mixed-layer depth observa- 
tions (not shown) exhibit the large-scale trends and much of the short- 
term variability shown in the simulation. 

4. Specification of atmospheric forcing 

The ocean mixed-layer model requires hourly estimates of wind speed, 
solar flux and total surface heat flux (if the salinity effects were in- 
cluded in the model, the precipitation rate would also be required). In the 
above examples, the atmospheric forcing was provided from 3 h observations 

at the ocean weather ships. It is thus an important question whether or 

not the atmospheric forcing can be provided at the required frequency over 

the ocean away from the weather ships (only a few of which remain in 

existence). A key premise in our research at the Naval Postgraduate 
School is that this atmospheric forcing can be derived from the Fleet 
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Numerical Oceanography Center (FNOC) atmospheric analysis/prediction 
models. Synoptic-scale wind fields are analyzed each 6 h and the 
atmospheric prediction models include hourly calculations of the surface 
heat budget over the ocean. Are these atmospheric model-derived fields 
adequate for ocean prediction? 

The present FNOC analyses of wind over the ocean make use of ship 
reports and satellite-derived cloud motion vectors. In many areas, there 
is adequate coverage to define the synoptic-scale variability. In other 
areas, the data coverage is poor, and there is reasonable doubt as to the 
validity of the analyses. A number of possibilities are being explored to 
enhance the observations of wind over the ocean. These include 
satellite-based instruments (such as on SEASAT) and over-the-horizon 
radar. If these instruments are to be useful for ocean prediction, they 
must be able to identify the oceanic regions of both high and low wind 
speeds. Given an accurate analysis of the wind field, the atmospheric 
model must provide the correct prediction of the winds with time. 

The components of the surface heat budget include the latent and 
sensible heat fluxes and the incoming and outgoing radiative fluxes. 

There are inadequate ship observations over the oceans to calculate these 
variables on the required (hourly) time scales. There are again a number 
of existing and proposed satellite instruments which may increase our 
capability to specify the heat flux components. One of the most difficult 
of the variables to observe remotely will be the latent heat flux, which 
depends on the near-surface specific humidity. For the present, we use 
the FNOC atmospheric model-derived heat fluxes (U.S. Naval Weather 
Service, 1975). The sensible and latent heat fluxes are calculated from 
bulk aerodynamic equations given the sea-surface temperature distribution. 
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Friehe and Pazan (1978) found good agreement between the FNOC winds and 
heat fluxes using independent observations over a two-week period. One of 
the more questionable aspects of the model heating package is the estimate 
of the cloud amount which is used in the radiation calculations. In the 
present formulation, cloud cover is simply related to the relative 
humidity in the column. 

An example of the time series of atmospheric forcing derived from the 
FNOC fields is shown in Fig. 6. The wind speeds are derived from the 6 h 
analyses using a cubic spline interpolation. Although the heat flux 
values are presently being archived at 6 h intervals, during the period 
shown the values were only available at 12 h intervals (details of the 
extraction and interpolation routines can be found in Gallacher, 1979). 

The wind and heat flux values derived from the FNOC fields appear to 
contain the synoptic-scale variability shown above to be important for 
ocean prediction. They do not contain the mesoscale variability that 
would be expected from point measurements. One measure of the suitability 
of these fields is the performance of the ocean model. 

5. Ocean thermal structure observations/predictions 

One of the primary reasons why ocean prediction has lagged behind 
atmospheric prediction has been the paucity of ocean observations. The 
global rawinsonde network with 12 h releases has been adequate, especially 
over land areas, to resolve the primary synoptic-scale atmospheric 
features. The time scales of the oceanic flow are considerably longer 
than in the atmosphere. However, the comparable space scales for 
baroclinic motion in the ocean are several orders of magnitude smaller 
than in the atmosphere. This requirement for observations on very small 
space scales will remain the greatest obstacle to ocean circulation 
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prediction. However, there are larger scale anomalous thermal structure 
features with space scales comparable to atmospheric cyclones. These 
ocean features are the object of research of the North Pacific Experiment 
(NORPAX). White and Bernstein (1979) have designed a ship-of-opportunity 
expendable bathythermograph (XBT) program called TRANSPAC for observing 
ocean thermal structure on space scales of a thousand or more kilometers. 
An example of the TRANSPAC data coverage for a month is given in Fig. 7. 
Although there are considerable gaps near the coasts, there are enough 
data in the central Pacific to define the thermal structure over a 15-20° 
latitude band. Haney (1980) has used the TRANSPAC data in simulation 
experiments with an ocean circulation model that are designed to test 
anomaly generation hypotheses. Because all TRANSPAC observations within a 
particular month are used in the analysis, this defines the time scale for 
the initialization and verification of the ocean model. 

A one-dimensional, mixed-layer model considers only the vertical 

fluxes of heat. Consequently a necessary condition for accurate model 

predictions is that the change in heat content from the ocean analyses 

must be nearly equal to the surface heat flux. A check (Elsberry, 

Gallacher and Garwood, 1979) of this condition using the TRANSPAC ocean 

temperature analyses and the FNOC surface heat fluxes is shown in Fig. 8. 

o 

The heat content in the upper 200 m along 175 W has been calculated 
relative to the 200 m temperature, which tends to remove the effect of 
vertically coherent fluctuations that may be related to non-mixing 
processes. Over most of the latitudinal band, there is a large decrease 
in heat content from September to November 1976. The vertical lines 
indicate the integrated total heat flux between 15 September and 15 
November 1976 calculated from the FNOC fields. North of 36°N there is 
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very good consistency between the two calculations. There is clearly 
disagreement between the surface fluxes and heat content changes south of 
36°N, but one can not determine which calculation is in error from these 
data alone. 

Simulations with the Garwood (1977) model in the two regions noted 

above are shown in Fig. 9. The initial (September) temperature profile at 

38°N, 175°W is rather unusual in that it has a mixed layer depth that is 

less than 10 m. The predicted October profile illustrates the large 

temperature decreases near the surface and temperature increases at the 

base of the layer expected with strong vertical mixing. Further cooling 

and deepening of the mixed layer is predicted from October to November and 

into December. The model prediction is in close agreement with the 

December TRANSPAC analysis. Agreement between surface fluxes and heat 

content change, plus the correct vertical temperature distribution, 

indicates that the parameterization of surface mixing processes is capable 

o 

of accounting for the thermal structure. A similar prediction at 32 N, 
175°W does not verify well. Although the largest discrepancy is found in 
the upper 100 m, the TRANSPAC values are consistently warmer than the 
model prediction throughout the upper 200 m. As discussed above in 
relation to Fig. 8, it is not clear whether the TRANSPAC analyses or the 
surface fluxes derived from the FN0C fields are in error in this region. 
There are two other possible explanations. The mixed-layer model may need 
to be adjusted for different conditions, or the one-dimensional models may 
not contain the necessary physics. In particular, a horizontal or 
vertical advective process that is not included in the model could 
possibly explain the deviations. 
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During the spring, the important feature to be predicted is the 
formation of the seasonal thermocline (see Figs. 4 and 5). The monthly 
TRANSPAC analyses have inadequate time resolution to verify the rapid 
transition predicted by the mixed-layer model. Examples of the mixed 
layer depth and temperatures at various longitudes along 38°N are shown in 
Fig. 10. These simulations with the Garwood model are driven by the time 
series of atmospheric forcing derived from the FNOC fields and are started 
from the March TRANSPAC analysis. The predicted mixed layer depth traces 
show considerable time and space variability. The transition to a shallow 
layer characteristic of summer conditions occurs before day 105 at 
longitude 215 (135°W). Apparent transitions occur at the other longitudes 
within a few days of this date - with some appearance of a lag in time as 
one proceeds eastward. Over the next 10 days, the mixed-layer temperature 
begins to increase. However, increased mixing occurs subsequently that is 
sufficient to eliminate the warm, shallow layer between 175°E and 165°W 
(longitudes 175 and 195 in the figure, respectively). The accumulated 
heat is then spread over greater than 100 m and the surface temperature 
decreases. Around day 130, the mixed layer again shallows in the western- 
most traces and then remains above 60 m for the remainder of the period. 
The associated mixed-layer temperature traces increase rapidly following 
the transition to shallow mixed-layer depths (see Fig. 10b), much as in 
Fig. 4b. One can see from the space between adjacent temperature traces 
that the time of transition as well as response to the later forcing 
events causes horizontal variability along the latitudinal section. 

Similar variability is predicted along other latitudinal and longitudinal 
sections (not shown). 
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It is of interest to explore what features of the atmospheric forcing 
are most relevant to these predictions of spring transition. Because the 
transition date varies at each position, the atmospheric forcing was com- 
posited for 10 days prior and 20 days following the transition date. The 
resulting composites of mixed-layer temperature depth for the six longitude 
traces in Fig. 10 are found in Fig. 11. The day-to-day variability in 
mixed-layer depth prior to transition contrasts markedly with the smaller 
mean value and variability following the transition. One can also note a 
distinct change in the rate of mixed-layer temperature increase with time 
across the transition date. The composites of the forcing variables are 
shown in Fig. 12. It is clear that the mean wind speed on the day of tran- 
sition is much lower than during the 10 days prior to this date. Although 
the mean wind speed increases slightly during the next few days, the aver- 
age wind speed over the 20 days following the transition remains smaller 
than before transition. Such an extended period with less wind mixing (re- 
call this term is proportional to the cube of the wind speed) is required 
for the warm stable layer to become well established near the surface. The 
variations in solar radiation derived from the FNOC fields do not seem to 
be very well correlated with the transition date. The overall trend is 
toward increasing values. However, the maximum values appear to be 
modulated by a synoptic-period variability which is larger than the differ- 
ence in peak values from before to after the transition date. Similar com- 
ments apply for the total heat flux trace, except that the nighttime 
upward heat fluxes are smaller after transition. This is consistent with 
the smaller wind speeds found during this period. One then finds a trend 
toward larger net downward heat flux across the transition date, which 
contributes to the increased warmer temperature. However, it appears from 
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these composites that the major factor in the transition is the distinct 
decrease in wind speed for an extended period of time. 

It should be emphasized that the net downward flux of heat is a neces- 
sary requirement for the occurrence of a transition, and the establishment 
of a seasonal thermocline. A cross-section of mixed-layer depths and tem- 
peratures along 175°W during the spring as predicted by the Garwood model 
is shown in Fig. 13. There is again a pronounced reduction in variability 
in mixed-layer depth between the early and later periods for all latitudes 
north of 36°N. At these latitudes, the mixed-layer temperatures increase 
rapidly following the transition dates as was the case in Fig. 10b. How- 
ever, the mixed-layer temperature does not increase at latitudes south of 
36°N. In fact, the temperature at 30°N continues to decrease slowly 
throughout the period. The evolution of the temperature profile at 30N, 
175W is given in Fig. 14a. The initial profile is taken from the March 1976 
TRANSPAC analysis. It is clear that the model did not predict the forma- 
tion of a seasonal thermocline as indicated in the verifying data from the 
June 1976 TRANSPAC analysis. As noted above, the prediction was for con- 
tinued cooling of the upper layer. One explanation for this trend was the 
absence of a net downward flux in the derived atmospheric forcing (recall 
that an excessive upward heat flux in these latitudes was also found during 

the fall period in Fig. 8). An experiment was run in which the upward heat 

2 

flux was reduced by 10 cal/cm during each hourly time step. The model 
then predicted the development of a seasonal thermocline which is quite 
reasonable in terms of the verifying data (see Fig. 14b). It thus appears 
that relatively small bias in the total heat flux can have an important 
impact on the correct prediction of the magnitude of the seasonal thermo- 
cline. If it is indeed a bias in the FNOC heat flux, which is the problem 
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in these latitudes, it may not have serious consequences for the atmos- 
pheric predictions since the atmospheric model is not run for extended 
periods without new data for correcting the temperature and moisture 
profiles . 

The capability of an ocean model to provide time series of realistic 
thermal structure profiles given the correct boundary conditions is very 
important. An example is the long-term integration of the Garwood model 
with the ocean weather ship data shown in Fig. 5. In many regions of the 
ocean, there will be long periods without new observations. Without a 
prediction model, the only option is to revert toward a climatological 
profile. If one had confidence in the calculations of the local heat 
budget and the wind forcing, it would be possible to use the mixed layer 
model to update continually the ocean thermal structure until new observa- 
tions obtained. It appears that an analysis-prediction-analysis cycle, as 
used for the atmosphere, would be a useful alternative to a system that 
simply reverts to climatology in the absence of new observations. This 
approach would only apply in oceanic regions where vertical mixing pro- 
cesses are dominant. In other regions, the prediction model would have to 
include advective effects. 

6. Justification for ocean prediction models 

Given that one has the potential to predict some phenomena, it is 
still necessary to justify economically the costs of developing and 
running the model. In the case of weather prediction, there are 
tremendous economic benefits because of the direct effects on the lives of 
people everywhere. The justification may be less dramatic in the ocean 
case, but a few possibilities should be mentioned, (i) Knowledge of the 
heat storage and distribution in the upper ocean may help us understand 
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our climate on time scales beyond a few weeks. Whereas we have only 
considered here the response in the upper ocean to atmospheric forcing, 
understanding our climate may involve feedback between the ocean and 
atmosphere, requiring coupled atmospheric-ocean models. Demonstrating 
that the ocean prediction model can correctly predict the ocean response 
to given atmospheric forcing, however, seems to be a necessary first step 
to interactive models. (ii) The National Marine Fisheries Service could 
use knowledge of the ocean thermal structure to improve fisheries manage- 
ment. (iii) The U.S. Navy is a primary customer for an ocean prediction 
model because the detection of enemy submarines by acoustical methods 
requires knowledge of the ocean thermal structure. The Navy has recently 
organized the Naval Oceanographic Research and Development Activity 
(NORDA) to provide the link between ocean research and the application in 
the fleet. The Numerical Modeling Group at NORDA has been tasked to test 
and develop an ocean analysis and prediction model. With the establish- 
ment of this group, and the anticipated computer upgrade at FNOC, which 
will be required to run such a model, it appears that the Navy has the 
organization structure and resources to bring an operational ocean 
prediction model into reality. 

7. Summary 

We have used the experience in numerical weather prediction as a 
framework for discussing the potential development of a limited ocean 
prediction capability during the 1980 f s. We have used some recent 
research results at the Naval Postgraduate School to illustrate some of 
the aspects that must be treated if an ocean model is to be integrated 
with real data and forcing. 
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(i) An ocean model that is to be driven with real atmospheric 
forcing must respond on a wide range of time scales. Ocean circulation 
models with embedded mixed layers are being developed by several groups 
for this purpose. The ability to treat both advective and mixing 
processes will be a prerequisite for predicting phenomena such as ocean 
fronts . 

(ii) A knowledge of the wind forcing during storms as well as during 
extended periods of low wind speeds is necessary to predict properly signi 
ficant changes in ocean thermal structure. We presently have analyses of 
the synoptic-scale variability in the wind forcing over the shipping lanes 
Further information from the data-sparse regions may be derived by remote 
sensing systems. In the predictive mode, the ocean model will be limited 
to the length of time that accurate wind fields can be predicted. 

(iii) Over large regions of the oceans there is an approximately 
local heat balance. This will permit "first-generation" models that are 
one-dimensional. Various ocean mixed-layer models have the capability of 
predicting the first-order changes in the ocean thermal structure given 
the correct forcing. These models need to be compared over a range of 
time scales and ocean conditions to select a candidate for the first- 
generation model. 

(iv) The mixed-layer models require a specification of the solar 
radiative as well as the total heat flux during the prediction period. It 
is proposed that the heating package of the atmospheric prediction model 
be used as an indirect representation of the thermal forcing. Further 
comparisons, such as in Friehe and Pazan (1978), of the model-derived heat 
flux components with point observations should be made. 
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(v) One benefit of an ocean prediction system lies in providing a 
better representation of the existing ocean thermal structure in data- 
sparse regions. An analysis-prediction-analysis cycle will carry forward 
information from limited observations. If the ocean model includes 
advective processes, the information from data-rich areas will be 
propagated into adjacent regions. 

It is anticipated that the development of an ocean prediction model 
will reveal many shortcomings in our knowledge of ocean processes. This 
will provide a stimulating environment for research and development 
efforts in all areas related to an ocean analysis and prediction model. 
Likewise, we will learn more about the data requirements necessary for 
accurate predictions and thus be able to deploy our limited resources more 
efficiently. If the experiences of numerical weather prediction are 
indicative, the decade of the 1980's should be an exciting period as ocean 
prediction models are developed. 
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Figure 1. Ocean temperature distribution predicted using an embedded 

mixed layer model within an ocean circulation model. Cross- 
section is aligned along the track of the model storm which 
was located at 450 km at the initial time and has moved toward 
the right at 5 m s -1 and is at 1319 km after 48 h. 
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PERCENTAGE OF OBSERVATIONS 



3 

Figure 2. Cumulative percentage of rank-ordered 3 h u#, u# and Q g 

values versus cumulative percentage of observations based on 
24-, 15- and 23-year samples (September to December) at OWS P, 
V and N, respectively. 
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Figure 3. Daily-averaged values of (A) sea-surface temperature, ( B) 

mixed layer depth, (C) insolation (Q ) and surface heat flux 

s 

*3 

plus back radiation (Q ), and (D) u, at OWS P during 1963. 

a 

Dashed lines represent the long-term mean (2*1 year) values. 
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Figure 4. Similar to Fig. 3 except at OWS V during 1959. 
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Figure 6. Wind speed, total heat flux and solar heat flux at 32N, 175W 

during 1977 extracted from FNOC atmospheric model analyses and 
12 h predictions. 
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Figure 8. Heat content (10 cal cm ) relative to the 200 m temperature 
calculated along 175°W from TRANSPAC analyses in September 
(solid) and November (dashed) 1976. Vertical bars indicate 
the cumulative surface heat flux between 15 September and 15 
November 1976. 
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Figure 9. Initial (circles) temperature profiles at (a) 38°N, 165°W and 
(b) 32°N, 175°W from September 1976 TRANSPAC analysis and mean 
predicted values at 10 m intervals for months of October 
(triangle), November (horizontal dash) and December 1976 
(cross). Verification data from December 1976 TRANSPAC 
analysis are given at 0, 30, 60, 90, 120, 150 and 200 m 
(diamond) . 
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Figure 10. Model-predicted layer (a) depth and (b) temperature during 
1976 along 38N. Initial spacing between the traces at each 
longitude are 100 m and 2 °, respectively. 
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Figure 11. Composite of model-predicted mixed layer temperature (upper) 
and depth (lower) for 6 longitudes as in Fig. 10, relative to 
spring transition date (zero day). 
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Composite of atmospheric forcing parameters relative to spring 
transition date as in Fig. 11. Wind speed, total heat flux 
and solar heat flux are extracted from FNOC atmospheric model 
analyses and 12— h predictions. 
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Figure 13. Model-predicted mixed layer depth and temperature as in Fig. 
10 except along 155W during 1976. 
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Figure 14. Model predictions as in Fig. 9 except at 30N, 175W and for 
March (circle), April (triangle), May (horizontal dash) and 
June (cross) 1977. Verification data are from June 1977 

TRANSPAC analysis, (a) Unmodified surface heat flux and (b) 

-2 -1 

reduction of upward heat flux by 10 cal cm h 
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